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Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  he  converted  to  SI  units  as  follows: 


Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

pounds 

0.45359 

kilograms 

Conversion  Factors,  Non-SI  to  SI  Units  of  Measurement 


vu 


1 


Introduction 


Description  of  Problem 


People  are  often  trapped  inside  collapsed  buildings  in  the  aftermath  of  a 
catastrophic  event  such  as  an  earthquake,  tornado,  or  explosion.  A  trapped  person’s 
ability  to  call  out  for  help  is  greatly  diminished  due  to  dehydration  and  weatoess  after 
several  days  of  being  trapped.  Many  victims  of  such  events  try  to  signal  to  rescuers 
by  hitting  or  tapping  on  the  ground  or  parts  of  the  structure.  One  of  the  techniques 
rescuers  use  to  detect  distress  calls  is  seismic  geophones  or  accelerometers.  These 
detect  the  elastic  waves  generated  by  the  trapped  person(s),  which  travel  through  the 
building  debris.  The  ability  of  rescuers  to  hear  these  "thumping"  sounds  (distress 
calls)  is  hindered  by  the  presence  of  background  noise.  These  background  noises 
include  mechanical  devices  such  as  bulldozers,  cranes,  or  other  traffic.  Urban 
enviromnents  also  produce  a  generally  high  level  of  coherent  background  "noise", 
generated  by  machinery,  vehicles,  etc.  which  are  in  continuous  operation.  The  goal 
of  this  effort  is  to  improve  the  quality  of  these  elastic  wave  signals  by  digitally 
processing  them  to  remove  the  background  noise. 

This  report  discusses  several  methods  of  separating  these  "thumping"  signals  from 
urban  "noise".  Both  simulated  and  field  data  are  presented  to  assess  the  potential 
strengths  and  limitations  of  each  method.  Several  techniques  of  hardware 
implementation  are  presented  and  recommendations  for  further  research  are  discussed. 


Description  of  System 


The  two  most  important  considerations  that  must  be  investigated  for  a  prospective 
noise  cancellation  system  are  portability  and  speed.  The  supporting  hardware  must  be 
easily  moved  and  operated  in  the  field.  Generally,  it  must  be  powered  with  batteries 
and  draw  little  current.  Any  noise  reduction  or  cancellation  signal  processing 
algorithms  must  operate  in  real-time.  A  typical  system  with  the  above  attributes 
would  have  the  following  characteristics: 

a.  It  should  be  less  than  1  cu  ft  in  size  and  weigh  less  that  15  lb. 
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b.  It  should  be  battery  operated  with  removeable  rechargeable  batteries  that 
power  the  system  for  at  least  6  hours. 

c.  It  should  digitize  the  outputs  of  the  geophones  and  process  the  data  in  real¬ 
time.  The  result  should  be  signal  conditioned  and  channeled  to  earphones  on 
which  rescuers  could  listen. 

Rescuers  need  to  arrive  at  a  disaster  site  and  deploy  the  listening  systems  at 
locations  where  people  are  believed  to  be  trapped.  Quick  determinations  need  to  be 
made  as  to  whether  people  were  communicating  distress  signals  by  voice  or  tapping. 
Generally,  there  is  not  time  for  detailed  digital  signal  processing,  spectral  plots,  or 
other  custom  signal  conditioning  in  order  to  evaluate  results.  Portability  and  real-time 
processing  are  imperative. 
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Methods  of  Addressing  Noise 
Cancellation 


Noise  cancellation  was  addressed  by  several  methods  of  digitally  processing 
previously  collected  data.  Finite  Impulse  Response  (FIR)  and  Infinite  Impulse 
Response  (DR)  filters  may  be  designed  to  band-pass  a  limited  bandwidth  of  data  to 
reduce  background  noise.  Noise  Cancellation  techniques  are  also  applicable  where  the 
developed  algorithm  automatically  designs  a  digital  filter  which  cancels  or  reduces 
background  noise.  Beamforming  has  also  been  applied  to  specifically  reduce  noise 
coming  from  all  but  one  direction. 

To  determine  the  best  method  of  signal  processing  to  use,  the  characteristics  of 
both  the  signal  and  the  noise  need  to  be  considered.  The  anticipated  thumping  signal 
is  a  set  of  narrow  band  signals  (short  impulses)  between  0.5  and  5.0  seconds  apart. 
These  signals  are  generated  by  a  person  hitting  a  section  of  a  collapsed  building  or 
ground  and  causing  it  to  vibrate  at  a  resonant  frequency.  The  background  noise  will 
have  a  more  varied  spectrum.  Some  machinery  may  produce  narrow  band  signals 
while  operating  at  a  constant  RPM,  but  as  the  speeds  and  loads  endured  by  these 
machines  change,  so  will  the  received  signal  spectrum  be  modified.  Other  noise 
sources  may  produce  signals  with  a  much  wider  bandwidth.  The  combinations  of 
these  noise  sources  generally  provide  a  very  complicated  noise  spectrum. 

Digital  Filtering 


Digital  filtering  provides  a  very  simple  method  of  processing  the  data.  If  the 
narrow  bandwidth  of  the  signal  is  known,  this  portion  of  the  spectra  can  be  bandpass 
filtered  and  a  great  portion  of  the  background  noise  will  be  eliminated.  This  is  the 
easiest  of  all  methods  considered  to  implement.  However,  for  proper  application  the 
bandwidth  and  center  frequency  of  the  signal  must  be  establish^.  This  information 
would  not  normally  be  known  in  advance.  Additionally  each  site  condition  has  its 
own  unique  characteristics  which  cause  the  bandwidth  and  center  frequency  to  differ 
from  locale  to  locale.  However,  one  practical  option  is  to  divide  the  spectrum  into 
multiple  narrow  bands  and  band  pass  each  section.  The  output  of  each  section  could 
be  monitored  in  order  to  detect  the  distress  "taps."  Certain  sections  of  the  audio 
spectrum  have  been  clearly  ruled  out  as  containing  survivor  response  taps.  These 
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may  be  always  filtered  out  of  the  collected  signals.  For  example  signals  under  20  Hz 
ne^  not  be  processed  since  they  carmot  be  heard  by  the  human  ear.  It  is  unlikely  that 
signals  above  250  Hz  would  originate  from  a  tapping  signal  source.  Consequently  a 
set  of  nine  or  more  filters  between  20  and  250  Hz  with  overlapping  pass  bands  would 
provide  excellent  separation  of  the  signal  from  background  noises  outside  the 
passband.  The  user  would  however  have  the  problem  of  monitoring  all  nine  bands  in 
order  to  determine  if  a  distress  call  was  present,  or  to  determine  which  band  contained 
the  greatest  signal  to  noise  ratio. 

A  modified  portable  computer  could  be  utilized  in  monitoring  all  nine  chaimels  and 
assisting  the  listening  system  operators.  This  would  be  accomplished  by  scarming 
each  band  as  it  is  filtered  and  numerically  searching  it  for  the  high  peaks  or  impulses 
produced  by  a  survivor’s  tapping.  The  determination  of  which  signals  constitute  an 
impulse  should  be  made  by  measuring  the  signal’s  peak  value  relative  to  its  overall 
noise  level.  The  algorithm  in  the  computer  could  then  select  the  frequency  band  that 
separates  the  data  and  noise  most  noticeably.  By  this  process  the  survivor’s  signal 
will  be  robustly  and  optimally  filtered  for  each  varying  site  condition  and 
environment. 

Either  HR  or  FIR  filters  could  be  implemented  for  use  in  such  systems.  FIR  filters 
offer  the  security  of  assured  numerical  stability,  but  HR  filters  can  be  designed  to 
operate  more  quickly.  These  trade-offs  must  be  carefully  considered  and  investigated 
when  designing  a  final  implementation  for  a  trapped  survivor  detection  and  location 
system.  Both  speed  and  stability  are  of  critical  importance  since  real-time  operation  is 
required. 

FIR  filter  with  the  transfer  function 


M 

H{z)  = 


x:=o 


(1) 


can  be  implemented  using  the  following  difference  equation: 


u 

y(n)  = 

i=0 

Where  y(n)  is  the  filtered  output,  x(n)  are  the  input  samples,  h(k)  are  the  filter 
weights,  and  M  is  the  number  of  filter  points. 
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IIR  filters  with  the  transfer  function 


Hiz)  = 


t=0 


N 


k=l 


(3) 


can  be  implemented  using  the  following  difference  equation: 


N  u 

)<«)  =  5^  a*  3<«-^)  +  E 

i=l  Jfc=0 

Where  "y(n)"  are  the  output  samples,  "x(n)"  are  the  input  samples,  and  "a"  and  "b" 
are  the  fflter  coefficients. 

Both  HR  and  FIR  filters  are  discussed  in  detail  in  Oppenheim  and  Schafer  [2]. 

Noise  Cancellation 


Noise  cancellation  can  also  be  employed  in  the  form  of  an  adaptive  fiilter. 
Conventional  filtering  accomplishes  noise  cancellation  in  an  open-loop  fashion  with 
FIR  or  HR  filters.  Adaptive  filters  operate  in  a  closed-loop  system  where  the  output 
is  monitored  to  determine  the  error  of  the  system.  This  error  is  then  used  to  modify 
or  change  the  adaptive  filter.  Several  different  algorithms  have  been  proposed  to 
implement  adaptive  filters,  however,  most  use  a  form  of  the  Least  Mean  Square 
(LMS)  algorithm  presented  by  Widrow  [1].  A  closed-loop  adaptive  filter  is  shown  in 
Figure  1.  In  this  example  an  input  signal,  s,  is  supplied  to  a  processor  and  an  output, 
y,  is  obtained.  This  output  is  then  compared  with  a  desired  output,  d,  and  the 
performance  or  error,  e,  is  measured.  This  error  is  then  used  by  an  adaptation 
algorithm  to  adjust  the  processor  until  a  minimum  amount  of  error  is  present.  As  the 
characteristics  of  the  input  and  desired  ouq)ut  change  so  will  the  learned  transfer 
function. 

If  there  is  a  linear  relationship  between  s  and  d,  the  adaptive  filter  will  learn  it  and 
the  output  e  will  show  the  signal  with  the  linearly  correlated  signal  removed. 

Single-channel  noise  cancellation 


A  single-channel  noise  cancellation  system  can  be  implemented  in  the  form  shown 
in  Figure  2.  In  this  diagram  a  signal,  s,  is  used  as  the  desired  signal  input.  A 
delayed  version  of  this  signal  is  then  used  as  the  main  input  to  the  adaptive  processor. 
If  this  input  signal  contains  periodic  signals,  both  the  initial  signal  and  its  delayed 
version  will  be  well  correlated.  The  adaptive  processor  will  then  learn  a  transfer 
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Figure  1 .  Closed-loop  adaptive  filter 


fuaction  that  predicts  what  future  values  of  the  signal  will  be.  The  adaptive 
processor’s  output  and  the  delayed  initial  signal  are  subtracted  to  determine  the  error. 
If  the  initial  signal  is  very  periodic  the  error  will  eventually  decrease  to  near  zero. 
However,  as  noncorrelated  data  occurs  in  the  signal  it  can  be  observed  at  the  error 
output,  e. 


This  form  of  noise  cancellation  can  be  used  in  applications  in  which  the  noise  is 
periodic  and  the  desired  signal  is  not.  One  example  of  such  a  situation  is  when 
human  voices  are  being  detected  in  the  presence  of  machine  noise  [8]. 
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Two-channel  noise  cancellation 


The  basic  form  of  a  two-channel  noise  cancellation  system  is  shown  in  Figure  3. 

In  this  diagram  one  input  is  a  desired  signal,  s,  combined  with  noise,  n,  to  form  s+n. 
The  other  input  contains  just  noise  that  is  correlated  with  the  noise  in  the  first 
channel.  The  signal  and  noise  channel  is  used  as  our  desired  output  and  the  correlated 
noise  channel  is  used  as  our  input  to  the  adaptive  process.  The  adaptive  process  will 
eventually  learn  the  linear  relationship  between  the  noise  in  the  first  channel  and  the 
noise  in  the  second  channel.  As  this  relationship  is  learned  the  output  of  the  adaptive 
process  will  become  equal  in  magnitude  and  phase  with  the  noise  in  the  first  channel. 
Since  they  are  subtracted  to  obtain  the  error,  this  difference  will  also  be  equal  to  the 
signal  without  the  corrupting  noise. 


Figure  3.  Two-channel  noise  cancellation 


Two-channel  noise  cancellation  systems  have  been  used  in  applications  such  as 
reducing  the  noise  in  the  cockpit  of  an  airplane,  or  in  detecting  Ae  ECG  of  a  fetus 
when  it  is  corrupted  by  noise  from  the  mother’s  heart  beat  [9]. 

The  LMS  Algorithm 


The  Least  Mean  Squared  (LMS)  algorithm  is  the  equation  used  by  the  one  and  two 
channel  noise  cancellation  systems  to  modify  the  filter  weights  and  thus  continually 
redraign  the  filter’s  response.  This  enables  it  to  adapt  to  a  changing  noise 
environment.  The  LMS  algorithm  is  derived  and  discussed  in  detail  in  Adaptive 
Signal  Processing,  by  Widrow  and  Steams  [1].  The  adaptive  filter  or  linear  combiner 
can  be  configured  as  a  transversal  filter  as  shown  in  Figure  4.  The  output  is  then 
calculated 
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(5) 


Y  = 


xl 


The  error  can  then  be  calculated 

=  d,  -  Xl 


(6) 


Then  the  filter’s  weights  can  be  adjusted  with  the  formula 

h;.,  (7) 


As  can  be  seen  from  formula  7,  the  weights  are  continually  being  adjusted  in 
proportion  to  the  error  e.  It  is  important  to  note  that  if  both  &e  desired  signal  and  the 
input  signal  are  the  same,  no  error  will  result  and  no  change  in  filter  weights  will 
occur.  Also,  if  there  is  not  a  linear  relationship  between  the  two  signals,  the  weights 
will  change  at  random  and  no  noise  reduction  will  occur. 
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Figure  4.  Linear  combiner  as  transversal  filter 
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3  Simulated  Data 


In  order  to  test  the  proposed  digital  processing  algorithms,  simulated  data  was 
generated.  It  was  used  with  the  proposed  algoridims  to  evaluate  how  well  it  worked 
under  controlled  conditions. 

Description  of  Data 


The  simulated  data  consisted  of  eight  tests.  The  simulations  used  four  noise 
sources  and  one  signal  source.  Two  noise  sources  were  correlated  and  two  were 
uncorrelated  noise.  A  diagram  of  the  model  for  the  simulation  is  shown  in  Figure  5. 

This  model  contains  a  wideband  noise  source,  Nl,  and  a  narrow  band  noise 
source,  N2.  Each  of  these  noise  sources  travels  through  a  different  transfer  function 
to  channels  1  and  2.  The  signal  also  travels  through  different  transfer  functions  to 
channels  1  and  2.  However  the  gain  of  the  transfer  hinction  between  the  signal 
source  and  channel  2  is  much  smaller  than  the  gain  between  the  signal  source  and 
channel  1.  In  addition  to  these  signals  random  noise  is  added  to  both  channels  1  and 
2.  The  transfer  fimctions,  Hl(z)  through  H6(z)  are  shown  in  Appendix  A.  For  each 
of  the  eight  tests  the  amplitudes  of  the  noncorrelated  random  noise,  the  amplitudes  of 
the  correlated  noise,  and  the  frequency  of  the  distress  signal  were  varied.  The  actual 
parameters  used  are  shown  in  Table  1. 

All  signals  were  normalized  to  a  1  Volt  RMS  amplitude  over  a  0.2  Sec  duration. 
All  wideband  signals  were  from  50  to  150  Hz.  All  narrow-band  signals  were  from 
58.9  to  61.1  Hz. 
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Table  1 

Simulated  Data  Parameters  of  a  1  Volt  RMS  Signal 

Test 

Signal 

Freq 

Hz 

Wide¬ 

band 

Amplitude 

Volts 

Narrow- 

band 

Amplitude 
58.9-61.1  Hz 
Volts 

Noise 

Amplitude 

Volts 

Test  SI 

111 

.1 

.2 

.01 

TestS2 

200 

.1 

.2 

.01 

Test  S3 

111 

.1 

.2 

.50 

Test  S4 

200 

.1 

.2 

.50 

Test  S5 

111 

.5 

.5 

.01 

Test  S6 

200 

.5 

.5 

.01 

Test  S7 

111 

.5 

.5 

.50 

TestS8 

200 

.5 

.5 

.50 

A  plot  of  each  channel  for  each  of  these  tests  and  its  spectrum  is  shown  in 
Appendix  A.  The  most  corrupted  test.  Test  S8,  and  the  least  corrupted  test.  Test  SI, 
are  shown  in  Figure  6. 
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Rgure  6.  Most  corrputed  (S8)  and  least  corrupted  (S1)  simulated  test 
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Filtering  Results 


To  test  the  ability  of  using  classical  FIR  filtering  to  detect  impulse  signals  a  set  of 
nine  FIR  band-pass  filters  was  implemented  using  61  points.  The  results  of  these  nine 
bands  for  test  S6  and  test  S7  are  shown  in  Figures  7  and  8.  As  can  be  seen  from 
these  plots,  the  filter  well  separates  the  impulse  signal  from  the  noise  so  the  bands 
containing  the  impulse  can  be  easily  identified.  The  nine  bandpass  FDR.  filters  were 
designed  with  the  following  bandwidths:  20-66  Hz,  43-89  Hz,  66-112  Hz,  89-135  Hz, 
112-158  Hz,  135-181  Hz,  158-204  Hz  181-227  Hz,  and  204-250  Hz. 

One  problem  with  this  method  is  that  the  user  would  need  either  advance 
knowledge  of  the  bandwidth  of  the  signal  to  be  detected  or  to  listen  to  all  nine  bands 
in  order  to  determine  which  band  has  the  best  signal  to  noise  ratio.  A  better  method 
of  accomplishing  this  is  to  allow  a  program  or  algorithm  in  a  modified  PC  computer 
to  determine  which  band  has  a  signal  that  most  resembles  the  distress  calls.  This  can 
be  achieved  by  calculating  the  ratio  of  the  maximum  value  in  each  band  to  the 
standard  deviation  of  the  continuous  noise  level  of  each  band.  This  will  give  a  near 
estimate  of  the  signal  to  noise  ratio  (SNR)  in  each  of  the  bands.  The  filtered  spectral 
sections  with  the  highest  SNR  can  then  be  chosen  for  continuous  monitoring  at  that 
specific  site. 

In  order  to  demonstrate  how  well  this  method  works,  each  of  the  eight  simulated 
tests  were  passed  through  a  set  of  bandpass  filters  and  the  one  with  the  best  SNR  was 
chosen  for  comparison  with  the  original  signal.  The  results  of  these  tests  are  shown 
in  Figures  9  through  12.  It  should  be  noted  when  viewing  these  plots  that,  it  is  not 
the  signal’s  amplitude  that  is  important  but  the  signal’s  amplitude  relative  to  the 
background  noise  level. 

The  results  show  that  considerable  noise  reduction  can  be  realized  by  this  type  of 
filtering.  In  the  even  numbered  simulated  test,  the  signal  and  wideband  noise  were  at 
distinct  different  frequencies,  but  in  the  odd  numbered  tests  the  signal  frequency  was 
near  the  center  frequency  of  the  wideband  noise.  It  can  be  seen  from  Figure  7  that 
when  the  signal  and  noise  are  separated  in  frequency  excellent  noise  reduction  is 
achieved.  Filter  8,  for  example,  shows  a  much  greater  SNR  than  the  raw  data  in  test 
S6.  It  should  also  be  noted  Aat  when  the  signal  and  noise  are  not  separated  in 
frequency,  less  impressive  noise  reduction  is  realized.  This  can  be  observed  in  Figure 
8  ^est  S7)  where  some  noise  reduction  occurs  when  comparing  Filter  4  with  the  raw 
data.  However,  the  noise  reduction  is  much  less  than  that  achieved  in  Test  S6. 
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Figure  7.  Bandpass  Filtering  Test  S6 
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Figure  9.  Rltered  Output  Test  SI  and  S2 
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Hgure  1 1 .  Filtered  Output  Test  S5  &  S6 
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This  demonstrates  that  for  the  bandpass  filter  to  work  well,  the  signal  and  noise  must 
be  separated  in  frequency. 
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Noise  Cancellation  Results 


Both  the  single-channel  (predictor)  and  the  two-channel  (noise  cancellation) 
routines  were  used  to  process  the  simulated  data.  For  the  single-channel  noise 
cancellation  routine  to  function,  the  input  must  be  very  periodic.  Channel  1  was 
delayed  by  50  msec  and  then  used  with  itself  to  filter  out  periodic  background  noise. 
The  results  for  the  eight  simulated  tests  are  shown  in  Figures  13  through  16. 

The  two-channel  noise  cancellation  routine  was  also  implemented  ^d  used  to 
process  the  eight  simulated  tests.  For  this  method  to  work  the  noise  in  both  channels 
must  be  correlated  and  the  signal  must  not  be  correlated.  The  results  are  shown  in 

Figures  17  thru  20.  _ 

The  results  show  that  both  single-channel  and  two-channel  noise  cancellation 
techniques  provide  good  noise  reduction  with  the  simulated  data.  For  the  two-channel 
noise  cancellation  to  work  effectively,  the  noise  on  the  two  channels  must  be 
correlated,  or  for  one-channel  noise  cancellation  to  be  effective,  the  noise  must  be 
periodic  (correlated  with  itself).  In  simulated  Tests  SI,  S2,  S5,  and  S6,  the  noise 
signals  were  well  correlated,  but  in  Tests  S3,  S4,  S7,  and  S8  large  amounts  of 
noncorrelated  random  noise  was  added.  The  results  of  the  analysis  show  that  little  or 
no  SNR  improvement  is  realized  in  the  tests  with  noncorrelated  noise.  However 
substantial  improvement  is  realized  in  the  others. 
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Figure  13.  Single-channel  Noise  Cancellation  Test  SI  &  S2 
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Figure  14.  Single-channel  Noise  Cancellation  Test  S3  &  S4 
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Figure  16.  Single-channel  Noise  Cancellation  Test  S7  &  S8 
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Rgure  17.  Two-channel  Noise  Cancellation  Test  SI  &  S2 
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Figure  18.  Two-channel  Noise  Cancellation  Test  S3  &  S4 
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Field  Test 


Data  from  a  field  test  was  used  to  evaluate  the  effectiveness  of  these  noise 
cancellation  algorithms  on  actual  data.  Field  test  data  were  selected  and  digitized  at 
1000  samples/sec.  These  data  were  processed  similarly  to  the  previously  discussed 
simulated  tests. 

Description  of  Field  Test 


Field  tests  were  conducted  at  the  Mississippi  University  for  Women  in  Columbus, 
Mississippi.  A  device  to  simulate  a  thumping  distress  call  was  constructed  and 
placed  in  a  multistory  office  building  which  was  then  collapsed.  The  building  was 
then  demolished  in  a  manner  to  simulate  an  earthquake  induced  failure.  Signals  from 
the  "thumpers"  were  then  recorded  on  digital  tape  using  geophones  at  various  selected 
locations  near  and  on  top  of  the  building  debris.  The  tests  using  the  most  responsive 
geophones  were  digitized  and  processed.  These  tests  used  two  transducers  placed  .5, 
1.,  2.,  4.,  and  6.  feet  apart.  The  placement  of  sensors  is  a  very  critical  consideration 
in  noise  cancellation  problems.  If  the  two  sensors  are  too  close,  both  signal  and  noise 
will  be  correlated.  If  they  are  too  far  apart  the  transfer  function,  H(z),  between  the 
two  systems  may  become  too  complex  and  nonlinear  and  therefore  too  difficult  to 
learn. 

In  order  to  study  the  character  of  these  tests  the  sections  of  the  time  history 
containing  the  signal  (thumping)  were  identified  and  separated  from  the  raw  data. 

The  noise  and  signal  could  then  be  studied  separately  or  together.  Plots  of  the  Power 
Spectrum  of  each  channel  and  the  Coherence  between  each  channel  were  calculated 
for  the  signal,  the  noise,  and  the  raw  data  of  signal  plus  noise.  These  plots  are  shown 
in  Appendix  B. 

Field  tests  were  made  with  three  main  configurations  of  sensors.  The  three 
diagrams  are  shown  in  Figure  21. 
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Digital  Filtering  Results 


The  data  from  the  field  test  was  also  used  with  the  band-pass  filtering  scheme 
discussed  in  Chapter  3.  These  results  also  show  very  good  noise  reduction. 
Examples  are  shown  in  Figure  22. 

It  can  be  seen  from  the  spectral  plots  in  Appendix  B  that  the  noise  is  separated 
from  the  signal  in  frequency.  The  signal  is  approximately  90  Hz  and  the  noise  is  a 
narrow  band  around  150  Hz.  This  would  indicate  that  the  selective  digital  filtering 
method  should  provide  excellent  results. 
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Figure  22.  RItered  Output  Test  26  &  36 
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Noise  Cancellation  Results 


Both  the  single-channel  predictor  and  the  two-chaimel  noise  cancellation  algorithms 
were  used  to  test  their  results  on  the  field  data.  To  get  good  results  from  the  two- 
channel  noise  cancellation  algorithm,  it  is  required  that  the  noise  signals  be  correlated 
and  that  the  signal  not  be  strongly  present  on  both  channels.  This  was  not  the  case 
with  this  field  data.  The  signal  was  highly  correlated  on  each  of  the  channels.  It  can 
be  seen  from  the  spectral  plots  in  Appendix  B  that  both  signal  and  noise  are  strongly 
correlated.  To  get  good  results  from  the  single-channel  noise  cancellation  algorithms, 
it  is  required  that  the  noise  signal  be  highly  periodic.  Again  this  was  not  the  case 
with  this  field  data.  As  can  be  seen  in  Figure  23  and  24,  the  noise  cancellation 
techniques  produce  poor  results  with  this  data. 
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Figure  24.  Two-channel  Noise  Cancellation  Test  26  &  36 
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5  Comparison  of  Results 


When  the  signal  is  correlated  on  both  channels  the  output  of  the  filter  can  be  saved 
instead  of  the  error.  This  provides  a  time  history  of  the  correlated  signals.  If  the 
noise  is  uncorrelated,  the  results  wUl  give  improved  SNR.  However,  if  both  signal 
and  noise  are  correlated  or  if  neither  is  correlated  no  SNR  improvement  can  be 
realized. 

The  three  methods  of  processing  the  data  were  compared  to  try  and  determine 
which  method  produces  die  best  results.  In  order  to  make  a  comparison  in  a  form 
that  is  representative  of  the  volume  levels  a  person  might  hear,  the  data  was  converted 
to  a  decibel  scale.  This  was  done  by  processing  the  data  through  an  RMS  filter  and 
converting  this  amplitude  into  decibels.  The  following  formulas  were  used. 


OUTPUT^,  = 


1  " 

-E^f 

NU  ‘ 


(8) 


and 


Decibels  =  20  LOGKXpUTPUTj^ 


(9) 


This  measure  of  decibels  should  not  be  confused  with  the  exact  audio  decibels  that 
might  be  calculated  for  a  microphone.  But  the  change  in  decibels  should  be 
comparable  to  that  of  a  microphone.  Decibel  plots  for  all  three  methods  were  made 
and  compared  with  that  of  the  raw  data.  These  plots  for  all  simulated  tests  and  field 
tests  are  shown  in  Appendix  C. 

Data  was  extracted  from  these  plots  and  tabulated  in  Tables  2  and  3.  These  tables 
show  a  comparison  of  the  three  methods  used  to  process  the  data.  They  indicate  that 
the  Selective  Filtering  method  works  well  with  all  of  the  data.  The  single-channel  and 
two-channel  noise  cancellation  routines  work  well  when  the  levels  of  uncorrelated 
noise  are  low.  In  field  tests  where  the  signal  and  noise  were  both  well  correlated  the 
noise  cancellation  routines  gave  no  improvement  in  performance. 

These  results  indicate  that  the  "best-method"  of  filtering  differs  from  site  to  site. 
One  algorithm  works  best  when  the  signal  and  noise  are  s^arated  in  fi:equency,  while 
another  works  well  when  the  signal  and  noise  are  not  both  correlated.  How  well  the 
signal  and  noise  are  separated  in  frequency  and  how  well  they  are  correlated  are 
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functions  of  the  site,  the  transducer  placement,  the  structure,  and  the  type  of  urban 
noise.  The  responders  have  little  control  over  these  conditions.  The  "best-method  is 
therefore  one  which  can  select  between  multiple  different  methods  of  noise  reduction. 
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Table  2 

SNR  Levels  Before  and  After  Processing 


Raw  Data  One  Chan  Two  Chan  Filter  Test 

Ch  1  Noise  Can  Noise  Can  db 

db  db  db 


13.2894 

j - 

13.6256 

15.6528 

1 

28.4677 

l«— . . 

26 

20.5432 

17.4282 

20.3388 

31.4466 

28 

16.7020 

16.0571 

20.2530 

28.6329 

29 

18.8568 

14.8417 

16.4318 

29.4628 

30 

16.5906 

15.1096 

16.0942 

27.6939 

31 

18.0985 

16.3148 

17.0262 

27.6135 

32 

18.9397 

18.1548 

19.9984 

29.4657 

33 

17.7485 

18.7080 

20.3775 

30.6381 

34 

27.9245 

24.8290 

26.5698 

28.0173 

35 

20.5917 

19.0165 

19.9874 

24.4057 

36 

20.9519 

15.4994 

20.9057 

31.1359 

37 

20.8629 

17.6962 

20.6946 

30.4611 

38 

17.8353 

15.0043 

17.4284 

21.0796 

39 

19.6213 

18.0286 

22.1001 

26.5954 

40 

18.8190 

17.4334 

12.6307 

6.5276 

41 

23.5342 

27.8079 

28.7579 

33.7639 

SI 

23.0693 

30.9092 

34.6914 

29.3173 

S2 

11.7452 

10.8272 

11.3803 

19.9742 

S3 

11.8002 

10.7441 

10.4280 

19.8004 

S4 

13.2299 

25.8739 

34.8918 

32.9168 

S5 

13.5019 

25.7274 

31.3524 

47.6927 

S6 

9.3654 

7.4533 

6.6000 

14.5801 

S7 

9.7874 

8.9595 

9.1001 

19.9054 

S8 
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Table  3 

Signal  to  Noise  Improvement 

One  Channel 

Two  Channel 

Filter 

Test 

Noise  Can. 

Noise  Can. 

Noise  Can. 

No 

db 

db 

db 

0.3362 

2.3633 

15.1782 

27 

-3.1150 

-0.2044 

10.9034 

28 

-0.6448 

3.5510 

11.9309 

29 

^.0151 

-2.4250 

10.6060 

30 

-1.4808 

-0.4962 

11.1033 

31 

-1.7849 

-1.0723 

9.5149 

32 

-0.7849 

1.0586 

10.5277 

33 

0.9595 

2.6290 

12.8895 

34 

-3.0954 

-1.3547 

0.0927 

35 

-1.5752 

-0.6043 

3.8140 

36 

-5.4525 

-0.0461 

10.1840 

37 

-3.1667 

-0.1683 

9.5981 

38 

-2.8310 

-0.4069 

3.2443 

39 

-1.5927 

2.4788 

6.9740 

40 

-1.3855 

-6.1882 

-12.2913 

41 

4.2737 

5.2237 

10.2297 

SI 

7.8399 

11.6221 

6.2480 

S2 

-0.9180 

-0.3648 

8.2289 

S3 

-1.0560 

-1.3721 

8.0002 

S4 

12.6440 

21.6618 

19.6868 

S5 

12.2254 

17.8504 

34.1907 

S6 

-1.9120 

-2.7653 

5.2146 

S7 

-0.8279 

-0.6861 

10.1179 

S8 
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Additional  Signal  Conditioning 


In  addition  to  filtering  the  data  other  methods  of  improving  the  quality  of  the  data 
should  be  employed.  It  is  recommended  that  Automatic  Gain  Control  (AGC)  be 
employed  in  Ae  data  conditioning  process  if  the  multi-bandpass  filtering  scheme  is 
employed.  The  AGC  should  be  used  to  make  sure  that  the  signal  is  always  at  a 
proper  level  for  listening.  Additionally,  a  method  of  replaying  data  should  be 
employed  so  if  a  signal  is  heard  the  responder  has  the  option  of  replaying  a  section  to 
hear  what  is  thought  to  be  distress  signals. 

Automatic  Gain  Control 


AGC  is  accomplished  by  periodically  measuring  the  average  level  of  a  signal  and 
continually  adjusting  the  gain  so  it  is  always  at  ^proximately  the  same  sound  level. 
Considerable  discretion  must  be  used  in  determining  exactly  which  AGC  algorithm  to 
use.  If  all  signals,  no  matter  how  small,  were  amplified  to  the  same  level,  very 
distant  signals  might  be  received  and  confused  with  signals  from  the  location  being 
investigated. 

Delayed  Playback  Schemes 


Delayed  playback  schemes  can  be  easily  employed  with  today’s  portable 
computers.  Inexpensive  Random  Access  Memory  (RAM)  and  hard  disk  storage  can 
be  used  to  store  substantial  amounts  of  data  for  later  listening. 

Visual  Aids 


In  addition  to  hearing  the  signal  it  is  always  helpful  to  have  visual  displays. 

Active  matrix  color  displays  offer  the  user  the  ability  to  view  plots  of  time  histories  of 
data  which  often  reveal  qualities  of  the  signal  that  are  difficult  to  hear. 


Chapter  6  Additional  Signal  Conditioning 


41 


Beamforming 


The  signal-to-noise  ratio  of  a  data  acquisition  system  can  be  improved  by 
employing  Beamforming  techniques.  Beamforming  is  a  method  of  using  combinations 
of  sensors  to  provide  improved  results.  A  detailed  discussion  of  Beamforming  is 
found  in  Haykin  [7].  When  the  direction  of  an  approaching  wave  and  its  speed  of 
propagation  are  known,  delays  between  samples  of  channels  can  be  inserted  in  order 
to  align  each  channel  as  if  it  were  parallel  to  the  approaching  wave.  When  this  is 
accomplished  all  signals  reach  the  array  of  sensors  at  the  same  time.  An  average  of 
the  sensors  is  then  taken.  Since  the  steered  wave  is  in  phase  and  other  signals  are  out 
of  phase,  the  desired  wave  will  be  amplified  while  other  noise  is  attenuated. 

Although  Beamforming  seems  theoretically  applicable,  it  has  practical  problems 
that  prohibit  its  use.  Exact  placement  of  transducers  must  be  made  for  beamforming 
to  work  effectively.  In  the  awkward  setting  of  the  rubble  pile  of  a  collapsed  structure, 
fast  deployment  of  a  beamformer  array  is  not  feasible.  However,  some  advantage 
could  be  realized  by  deploying  additional  geophones,  accelerometers,  or  microphones 
so  extra  channels  could  be  processed  and  tested  for  distress  taps. 
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Potential  Hardware 
Implementations 


Several  methods  of  implementing  a  system  to  acquire  seismic  data,  filter  it  and 
reproduce  the  filtered  data  are  presented.  The  first  method  is  to  design  a  system  that 
is  unique  for  this  application.  Developmental  cost  for  this  method  would  be  high. 
However,  if  many  units  were  to  be  produced  the  total  implementation  cost  would  be 
low.  Another  method  would  be  to  equip  a  portable  programmable  computer  (PC) 
with  commercially  available  expansion  printed  circuit  cards  to  acquire,  process,  and 
ouq)ut  the  analog  signal.  Such  a  system  would  be  much  less  expensive  to  develop. 

Integrated  Circuits 


Rq)orts  of  many  systems  have  been  published  in  which  special  digital  signal 
processing  (DSP)  integrated  circuits  have  been  used  to  implement  digital  FIR,  DR,  or 
adaptive  filters.  Examples  of  these  integrated  circuits  include  the  Texas  Instruments 
TMS32010,  TMS3020,  TMS320C25,  the  Motorola  DSP56000,  the  Analog  Devices 
ADSP-2100,  and  the  NEC’s  uPD77230  [5].  These  chips  can  be  used  with  analog  to 
digital  conversion  (ADC)  integrated  circuits,  and  digital  to  analog  conversion  (DAC) 
integrated  circuits  to  acquire,  process,  and  r^roduce  a  filtered  signal  in  real-time. 
These  DSP  integrated  circuits  can  be  preprogrammed  to  perform  many  different 
algorithms. 

Printed  Circuit  Cards 


An  easy  system  to  implement  would  be  to  equip  a  portable  IBM  compatable  PC 
with  a  printed  circuit  card  that  performs  the  anaJog  to  digital  and  digital  to  analog 
conversions.  These  cards  are  available  from  many  different  vendors.  If  high  speed 
80386  or  80486  microprocessors  were  used  in  the  PC  the  computer  itself  might  be 
able  to  process  all  the  required  filtering.  However,  if  the  spe^  of  the  computer  is 
not  sufficient  to  perform  all  the  necessary  processing,  DSP  cards  are  currently 
avaUable  that  utilize  DSP  integrated  circuits  like  the  ones  discussed  previously. 
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Example  of  Proposed  System 


It  is  reconunended  that  the  noise  reduction  system  outlined  in  the  previous  section 
be  implemented.  This  system  is  being  proposed  as  a  prototype  and  developmental 
system  to  be  used  in  bo  A  simulated  and  actual  rescue  attempts.  Hardware,  that 
utilizes  PC’s  to  control  the  filtering,  offers  the  advantage  of  being  easily 
programmable.  It  also  has  the  advantage  of  providing  a  large  number  of  filtering 
options  to  the  user.  It  provides  storage  for  data  samples  that  may  be  acquired  and 
archived  for  future  study.  In  addition  to  the  real-time  operation  of  the  system,  it 
would  be  beneficial  to  save  an  occasional  sample  of  data  so  other  filtering  routines 
could  be  evaluated  on  real  data  after  the  emergency  response  has  ended. 

Once  the  routines  that  provide  the  best  results  are  determined,  the  practicality  of 
designing  an  even  less  expensive  system  utilizing  DSP  integrated  circuits  should  be 
considered.  A  typical  system  could  be  configured  as  follows: 

a.  Portable  IBM  computer  with  2  full  length  expansion  slots;  16  Mbytes  of 
RAM;  200  Mbytes  of  disk  storage;  Floppy  Disk  storage;  486  Processor 
running  at  66MHz;  Active  Matrix  Color  VGA  display; 

b.  National  instruments  AT-DSP2200  Dynamic  Signal  Acquisition  and  DSP 
Board  for  PC  AT.  [6] 

c.  Current  Corps  of  Engineers  STOLS  unit.  [10] 

The  DSP2200  can  provide  both  ADC,  DAC,  and  DSP  capabilities  in  one  card. 

All  this  equipment  can  be  commercially  purchased  with  maintenance  available  from 
the  vendor.  No  hardware  development  time  would  be  required,  and  the  software 
development  time  could  be  performed  in  a  high  level  language  like  FORTRAN  or 
PASCAL.  As  the  system  is  used  in  real  or  simulated  rescue  attempts,  the 
effectiveness  of  the  different  filters  could  be  evaluated  and  modified  as  necessary. 

The  DSP2200  operates  at  25  million  floating  point  operations  per  second 
(MFLOPS).  A  486DX  running  at  33MHz  can  perform  at  approximately  2.4 
MFLOPS.  One  FIR  filter  using  128  points  would  require  128  multiplications  and  128 
additions  or  256  floating  point  operations  for  each  sample  of  data.  If  data  was 
sampled  at  1000  samples  per  second  per  channel,  it  would  require  256,000  floating 
point  operations  for  one  channel  or  512,000  floating  point  operations  for  both 
channels.  This  indicates  that  at  least  4  different  filters  could  be  implemented 
simultaneously  on  the  486  computer,  requiring  4*(.512  MFLOPS)  or  2.048  MFLOPS. 
The  sampling  rate  of  1000  samples  per  second  and  the  filter  size  of  128  points  is  the 
maximum  that  would  be  required.  If  the  number  of  points  was  dropped  from  128  to 
32  the  total  number  of  FIR  bandpass  filters  could  be  quadrupled.  TTie  use  of  the 
DSP220  would  then  increase  the  complexity  or  number  of  filters  by  a  factor  of  10. 

This  proposed  system  will  provide  an  important  tool  for  evaluating  many  different 
noise  reduction  algorithms  in  both  simulated  and  actual  rescue  attempts.  It  should 
provide  considerable  improvements  over  existing  commercially  available  noise 
reduction  systems.  It  also  provides  considerable  flexibility  and  can  be  modified  as 
continued  research  reveals  more  information  concerning  noise  cancellation,  the  nature 
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of  the  distress  calls,  and  the  nature  of  the  site  environments. 
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8 


Conclusions  and 
Recommendations 


Three  ways  of  adaptively  filtering  data  have  been  presented  and  tested.  The  first, 
single  channel  adaptive  filtering,  assumes  that  the  noise  is  periodic  and  the  desired 
signal  is  not.  When  using  this  method  geophones  should  be  placed  as  close  to  the 
signal  source  as  possible.  The  second  method,  two  channel  adaptive  filtering, 
assumes  that  the  signal  received  on  one  channel  is  not  correlated  with  the  signal 
received  on  the  other  channel,  and  that  the  noise  received  on  one  channel  is  correlated 
with  the  noise  received  on  the  other.  When  using  this  method  one  geophone  should 
be  placed  near  the  signal  source  and  the  other  in  a  location  that  is  distant  from  the 
signal  source  but  in  the  same  general  location  as  the  noise.  The  third  method  assumes 
that  the  frequency  of  the  signal  is  different  for  that  of  the  noise,  and  that  the  character 
of  the  signal  is  known  to  be  an  impulse.  The  data  can  then  be  filtered  into  several 
frequency  bands  and  the  ratio  of  peak  signal  to  standard  deviation  will  provide  an 
estimate  of  how  well  the  noise  has  been  removed.  The  frequency  band  in  which  this 
ratio  is  highest  should  be  mrmitored.  Any  number  of  channels  can  be  used  with  this 
method. 

Results  show  that  the  third  method  provided  the  best  results  in  most  simulated  and 
field  tests.  In  the  field  tests  both  data  and  noise  were  correlated  so  well  that  the  two 
channel  adaptive  filtering  provided  no  improvement  and  usually  worsened  the  quality 
of  the  data. 

It  is  recommended  that  further  study  focus  and  include  the  implementation  of  a 
system  like  that  described  in  Section  7.3.  This  method  demonstrates  the  best  noise 
reduction  data  processing  method  as  applied  to  the  detection  of  trapped  victims. 
Preliminary  software  should  be  developed  to  provide  a  wide  variety  of  filtering 
options  as  addressed  in  this  study  and  be  adaptable  to  include  any  new  data  processing 
methodologies.  Data  should  be  obtained  from  additional  field  tests  and  processed  in 
real-time  to  demonstrate  the  robustness  of  this  signal  conditioning.  The  results  of 
each  of  the  filter  types  or  level  of  complexity  should  be  evaluated  for  its  effectiveness. 
Such  a  system  would  operate  substantially  better  than  the  current  commercially 
available  noise  reduction  systems.  It  should  be  deployed  to  field  crews  and  serious 
consideration  should  be  given  to  manufacturing  a  specialized  system  like  that 
described  in  Chapter  7. 
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Appendix  A 

Analysis  of  Simulated  Data 


Appendix  A  contains  plots  of  the  simulated  ground  transfer  functions  used  to  generate 
two  channels  of  simulated  data.  (Referenced  on  page  10). 

This  appendix  also  contains  the  simulated  data’s  time  history,  PSD,  coherence,  and 
transfer  function.  (Referenced  on  page  12). 
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Appendix  B 
Analysis  of  Field  Data 


Appendix  B  contains  plots  of  the  raw  time  history  for  two  channels  of  the  field  data. 
It  also  contains  the  raw  data’s  PSD,  coherence,  and  transfer  functions.  (Referenced 
on  page  32). 

This  appendix  also  contains  these  same  plots  for  its  signal  component  and  its  noise 
component.  (Referenced  on  pages  33  and  35). 
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Appendix  C 

Decibel  Plots  of  Simulated  and 
Field  Data 


Appendix  C  shows  plots  of  the  processed  results  for  both  the  field  test  and  the 
simulated  test.  The  plots  show  the  results  for  band-pass  filtering,  one-channel  noise 
cancellation  filtering,  two-channel  noise  cancellation  filtering,  and  no  filtering. 
(Referenced  on  page  38). 
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